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Abstract 

In [| the authors suggested a new algorithm for the numerical approximation of a 
BSDE by merging the cubature method with the first order discretization developed by 
Bouchard and Touzi Q| and Zhang [T(|. Though the algorithm presented In [5] compared 
satisfactorily with other methods It lacked the higher order nature of the cubature method 
clue to the use of the low order discretization. In this paper we Introduce a second order 
discretization of the BSDE in the spirit of higher order Implicit-explicit schemes for 
forward SDEs and predictor corrector methods. 
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1 Introduction 

The present paper Is concerned with the problem of numerical approximation to forward 
backward SDEs (FBSDEs henceforth). Let (CI, J 7 , P) be a probability space on which we 
have defined a triple of processes (X, Y,Z) which solve the decoupled forward-backward 
system: 

rt d t 

X t = X + V (X u )du + ^2 V i i X ") ° dW u 

,=1 t 6 [0, T] . (1 .1 ) 

Y t = <t>(X T ) + J f(X u ,Y u ,Z u )du- J Z u -dW u 

where W t is a d-dimensional Brownian motion and 

V k : R f ' ->■ R c ', k = d, f:l'xlx R d -> R 

are some appropriate functions. The system is called decoupled as the (backward) pro- 
cesses {Y,Z) do not appear in the dynamics of the forward component X. Systems of the 
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form (1 .1 I have received a lot of attention over the past twenty years primarily clue to their 
applications in the field of Mathematical Finance (see for example [6| and the references 
therein for details). 

Of equal Importance Is the fact that, the stochastic flow associated with I.e., the triple 
of processes (X^*', yM.Z^), (t,x) G [0, I]xi satisfying 

K* = X + f 5 V IK'") dLI f V ' {<'") dW u 

i=1 se[t,T], (1.2) 

V^ = <D(Xf) + J {(X^.Y^.Z^du- J Z!;X-clW u 

provides a Feynman-Kac representation for the (viscosity) solution of a class of seml-llnear 
PDEs. In particular, let u(t, x) : [0, T] x R d R d , be the viscosity solution of 

+ V(x) • Vu + ^Tr[V(x) V*(x)D 2 u] + f (t.x.u.VuV (x)) = 
u(T,x) = <t>(x) 

where V (x) = V (x) - 1^ =1 V V } (x)V s (x). In their seminal work Parcloux and Peng 
showed that 

u(t,x) = V/' x , (t,x) G [0, 7"] x R rf . (1.4) 

In addition, Ma and Zhang [14] showed that when this viscosity solution Is continuously 
cllfferentlable In Its spatial variables, we have the following representation for the solution 
of the PIPE 11.31 and Its gradient 



u(t,x) = K/' x , Z' t ' x = Vu(t,x)V{x), a.s. (1.5) 

From the perspective of numerical analysis, the above means that any probabilistic method 
for the resolution of (1.1) . provides an algorithm for the resolution of seml-llnear PDEs. As 
a result, the Interest In robust algorithms for their resolution Is high. 

Algorithms designed to solve the problem of numerical approximation of the backward part 
of {1 .1 1 consist of two parts: Firstly, the backward equation Is cllscretlzecl. This step Involves 
the use of one or more conditional expectations. Secondly, a numerical method Is grafted 
onto the chosen discretization to compute the conditional expectations Involved. Following 
the work of Bouchard and Touzl [3] and Zhang [ToT ]. the most popular discretization scheme 
Is given byQ 

Zf t( := J-E [ Y? tM AWi +1 \T tl ] , i = n - 1 (1 .6) 

0['+1 

Kft := E [ Yft w \T ti ] + f(t h X, , Yfa , Zft ) , i = n - 1 , 

where tt Is a given partition tt := {0 = fo < U < ■ • • < t„ = T} of [0, 7] and AW; + i = 
M/ t . +i - W ti , i = n-1. 



^ee Bender and Zhang 0] for an alternative approach based on Picard Iterations 
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Assuming that all coefficients of (1.11 are at least Lipschitz continuous in their spatial 
variables, we have, following [3J] and [1 61 ], that 



0<t<T L J L Jo 




(1.7) 



where {(Y",Zf), t > 0} are the step processes 



n-1 



(7-1 



Yt ■= Y. + ^U, W Zf := J^Z^ 1 [f ,,, +l) (f) + 1 t=t „. (1.8) 



In other words, the above discretization of the backward part achieves a convergence of 
order 1/n 1/2 , when n points are used, i.e. the same order as the strong convergence order 
of the Euler scheme for a (classical) SIDE. When more smoothness on the coefficients is 
assumed, it is shown in Gobet and Labart [8] that the rate of convergence of the processes 
(1 ,8| is of order 1/n. In fact an error expansion is obtained in j{| and the leading order 
coefficients in the error expansion are identified. 

To be more precise, the above results are proved for the case when the process X in (1 .61 is 
replaced by its Euler approximation X n and, respectively, the filtration { JF f; }" =0 is replaced 
by the natural filtration {J-*}" =0 associated with X". However the same proofs apply both 
to X n and X. 

To obtain an Lmplementable scheme using the discretization (1 ,6| . one has to develop a 
method for approximating the conditional expectation^] EfV^ \J~t] and ElYT^ AH//+1 \J- t ]. 
Various such methods have been introduced, based on Malliavin calculus [3], on projection 
on function basis [9], [10] and on quantization [1|. In Q, the authors suggested the application 
of the cubature method of [ijj], which is based on the ideas of Kusuoka jl2| . The overall 
rate of convergence of this second approximation step is still of order 1 fn " or of order 1/n 
when coefficients are smooth. 

In this paper, we introduce a discretization for the backward component of the BSDEs of 
order 1/n 2 . As for the discretization |1 ,6| , this new scheme will require the computation of 
conditional expectations. This can be clone by means of an order 5 cubature formula. Since 
the order 5 cubature method has local order error 1 In 3 , when combined with a second order 
discretization for BSDEs, will generate a genuine second order algorithm for BSDEs. 

To understand which terms such a scheme should include we revisit II .61 from the point 
of view of Stratonovich-Taylor expansions. Moreover, to deduce the rate of convergence of 
the method, we shall rely on the Feynman Kac representation 11 .5} and on some estimates 
on the bounds for the derivatives of (1 .3} obtained by Delarue [7|. 



2 Preliminaries 

Throughout the paper, we will use the following assumptions: 

2 When the process X In 11.61 is replaced by its Euler approximation, then EfV^f.^ and E[V 1 Ir t +l AH/; + i IJ 7 ^] 
need to be computed. 
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(A) The coefficients of the forward SDE K : M d -> R cl , i = 0,1 cl have all entries 

belonging to C^°(M. d ), the space of bounded infinitely clifferentiable functions with 
all partial derivatives bounded. We also assume that the matrix V := (V\\ . . . is 
elliptic. 

(B) The driver of the BSDE f:[0,1]xl l 'xMxM < '^l belongs to c]" ll2] ' m . The exact value 

for the parameter m shall be made precise as we proceed. 

(C1) The terminal condition <t> : K — * M. is Lipschitz continuous. 

(C2) The terminal condition <t> belongs to C," 1 (R C ';R). Again the value of m shall be deter- 
mined further on. 

We shall denote by K the constant that bounds all derivatives that appear in our assump- 
tions. When (C1) Is in force we shall assume that <t> Is /C-Lipschitz. Of course under (A), 
(B) with m > 1 and (C1) the system {1.1} has a unique solution such that 



To abbreviate notation we shall denote by Mi = M ti , t = 0, . . . , n, where M can be any of 
the processes that appear In this paper. Moreover we write E s [-] for E[-|.F S ]. Note that 
In the current (Markovlan) set up E[-|J r s ] = E[ |X S ]. We shall also consider as given a 
partition n :— {0 = to < fi < . . . < t n = T} of [0, T] and consider the related notation 

5 £ := t t - — ff-T , [^Wi'-Wi-W^, AWl := W! -Wl_ v / = 1 d, t = 1 n. 

Lastly, the driver of the BSDE shall be abbreviated as f(s,x) = f(x, u(s, x), Vu V(s, x)), 
where u is the solution of the semilinear PDE JT3J. 

Remark 2.1. Any attempt to discretize the backward part of |1 .1 ) should start with a 
method to produce (an approximation for ) the value of the forward diffusion at the times 
on the partition n. In this paper we assume that we have at our disposal the forward 
diffusion X at the times of the partition X t] , . . . ,X tn . This entails no loss of generality, 
as we aim to combine the second order discretization with the cubature method, which 
approximates directly the law of X at the times of the partition. 

For the benefit of methods that would combine the present discretization with a Monte 
Carlo simulation, that most probably requires some sort of discretization for the forward 
process X, we note that all results regarding the second order scheme (particularly Theorem 
\3.3\ and Corollaru \3.4i remain valid, once one replaces X with a second order approximation. 
A wealth of such approximation^ are presented on chapters 12- 15 offTfy. 

Working towards a higher order discretization of the backward part of (j .1} we shall rely 
heavily on the Stratonovich-Taylor expansions. Hence we need to fix notation and present 
some elementary facts regarding these expansions: 

3 Such approximations can be chosen to be explicit, in other words completely derivative free. 
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2.1 Stratonovich-Taylor expansions 



In the following we will use the set of multi-indices A = {0} UlJ^L^O, 1, c/}'" endowed 
with the norms | ■ | and ||-|| given by 

|/3| = length of (3, ||jS|| := |/3| + Carcllj : jSy = 0, 1 < 7 < |/S|}. 

Clearly |0| = |]0|] = 0. For a /3 = (j\,...,jt) e i we also write /3— = (ji . • • - »7/-i) and 
— /3 = (_/2i ■ • ■ .7;)- We define the subsets of A, 

A m = {/3eA: ||/3|| < m} and A' m = {/3 G A\{0, (0)} : ||j8|| < /77}. 

Given two multi-indices a = (ai , . . . , a^) and /3 = (/3-| , . . . , /3() we define their concatenation 
as a * /3 = (ai , . . . , a/t, j8i, . . . , j8/). For a suitably chosen function f and a multi-index 
/3 = (j8i, . . . ,/3;), we define the iterated Stratonovich/Lebesgue integrals as follows 

'f(s) |/3|=0 
W\t.s-= \ J* k-[f]t,udu l>,, j t = 0. 

f t s Jp-[f}t,u°clW'>(u) /> 1,7,^0 

Let U , j = 0, . . . , d be the operators: 

q q 
L> :=^Vfd Xk , y = 1 d, L° :=d t + ^2v k d Xk . 

fc=i /<=1 

The iteration of this family of operators is understood as : L a f := L" 1 . . . L°" f, a = 
(cr-i, . . . , a n ) and we also use the convention L?f = f. 

Given a multi index a = (a-[,...,a n ) the following identity is proven in Proposition 5.2.10 

of [El 

17 + 1 

W{j a [\\a, t = ^Jo.t (a1 a "^^K 7 = 1 d. (2.1) 

A direct consequence of (2J) (Corollary 5.2.11 of [ill) is that if a = (k, k k), k = 0, . . . , d 

with | or | = m, then 

JV}o,t = ^{J {k) W ,t) k (2.2) 

The (conditional) expectations of iterated integrals can also be computed. In particular we 
have a characterization for those iterated integrals that have non zero expectation (see, for 
example, j3j): 

Lemma 2.2. Let a = (it i r ) be an arbitrary multi-index with \\a\\ = m and t G [0, T]. If 

m is odd, then E[y a [1]o,f] = and if in is even then 



2 r -"'' 2 (m/2)! LI C G •A-m.r 

otherwise 



where A m r C A m is the set of multi-indices with a = ai * . . . * a,„/2 G A m , such that a; = (0) 
or a,- = (7,7), 7 G {1 d}. 
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A non empty subset Q C A Is called a hierarchical set If sup ae g \a\ < oo and —a E 
5 for any a £ 5\{0}. Given a hierarchical set Q we define the remainder set 

B(G) = {/3 e A\G : -PEG}. 

Note that A m and .A^ are hierarchical sets. If f : R" — > R Is a smooth function with all 
partial derivatives bounded, then by repeatedly applying the Ito-Stratonovlch formula, we 
obtain the Stratonovlch-Taylor expansion 

f(X^) = ^L a f(x)J a [\] , t + nL a f(-)]o,t 
=:^L a f(x)J a ['\} , t + R m (t,x,f) 

aeK 



for any hierarchical set /C. The second term on the right hand side of |2.3| Is called the 
m-th order remainder process and Ls denoted by R m (t,x, f). Obviously as in Increases, and 
for small t, the remainder process gets smaller and smaller. 

A hierarchical set that facilitates computations for Stratonovlch-Taylor expansions Is A(m). 
For this set we have that 

B(A(m)) = {(j)*P\j = d, ||jS|| = m}. 

To estimate the remainder that corresponds to B(A(m)), we need to estimate Iterated 
Stratonovlch Integrals of the form J a [g ( •, X . )] s f for appropriate function g such that the 
previous Integral makes sense. For < s < t with t — s < 1, we have 

sup E s \j a [g(-,X.)} SI ] < C(f-s)HI /2 ( sup {E s [\L a "g(r,X r )\ 1„^ ] 

s<r<t L ' J *s<r<t |2 4) 

+ E s [|g(r,X,-)|1 a „ =0 ]}) 

see chapter 5 of Kloeclen and Platen [TlJ. Applying the estimate |2.4| to every Index In 
the set B(A(m)) , provides us with the following estimate for the remainder of the Taylor 
formula, 

E[\R m (t,X,f)\ p } Vp < C(f-s) (m+1)/2 max sup E[ | L y f{r, X r ) \ p ] Vp . (2.5) 

»i<l|y||<'"+2o<r<f 



To conclude our preliminary results, we shortly discuss the backward PDE j1 ,3| . For the 
latter, when all coefficients are smooth and bounded, It Is well understood that It posseses 
a smooth and bounded solution. Moreover, a smooth solution u(t, x) will still exist when 
working only under (C1) for t £ [0, 7~). The solution Itself will be bounded on all [0, T]. Its 
derivatives however shall not. The next result characterises the behaviour of the derivatives. 

Theorem 2.3 (Delarue [7]). Let (A) and (13) hold true and assume further that <t> is Lips- 
chitz continuous. Then there exists a unique u £ 

c \ml2lm | ^ f) x Rd | thaf so[ves jj^j 

Moreover, for ang a £ A(m) we have 

||D» u (t,.)L<C (r ^l- 1)/2 , te[0,T), 
for a constant C independent of u. 
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In other words, when the final condition is Lipschitz continuous (and the driver is smooth), 
the derivatives of the solution of {1 ,3} , explode as t f T with the given rate. The exact rates 
are important to quantify the error of our method. 



3 One step discretization 



To understand the terms that a second order scheme should incorporate, we present first the 
intuitive arguments that lead to the discretization scheme |1 ,6| . We consider the backward 
part of between two successive times of the partition n 



f(X s , Y s ,Z s )ds 



f 



Z, ■ dW, 



and cliscretize the Riemann integral using the left hand side point, as in pj], thus leading 
to an implicit equation for Y t . and the stochastic part in the usual way, to obtain 



Y ti ~ Y tM + 5 M f(X ttl Y til Z ti ) - Z u ■ /\Wi 



+i • 



(3.1) 

By conditioning |3.1 | with respect to T tl in {3.1} we obtain a first order approximation for 

Yt, 



Y ti ~ E 



[V ti+1 |^ ti ] +5 M f(X tl ,Y tl ,Z tl ] 



(3.2) 



but for the presence of Zt r To treat the Z ti , we can m ultiply both sides of (33) by AH/' +1 , / = 
1, . . . , cl and condition with respect to Tt v to obtain 



Z 



E 



Y AM/ / +1 



5t 



+1 



/ = 1 d. 



(3.3) 



Observe that scheme j1 ,6| is just the backward iteration of equations |3.2| and {3.3 



Clearly, to achieve a higher order approximation of the BSDE, we need to add more terms 
in its (stochastic) expansion ideally without involving any derivatives. As far as the driver 
is concerned, the next elementary result tells us that the Crank-Nickolson rule indeed 
achieves a third order local error, thus leading to a second order global error: 

Lemma 3.1. Let assumptions (A) and (B) hold true, <t> £ Cu p ( K d ) and denote by f(t,x) = 
f(x,u(t,x),'Vu(t,x)V(x)), where u(t,x) is the classical solution to PDE {1 .3} . Then there 
exists a constant C independent of the partition and of u such that 



E, 



jT ' +1 f(s, X s )cls - ^±1 ( f(ti,X ti ) + f(t M , X tM ) ) 



max 

||a||=4,5 



,a u(ti 



+ 1. 



Proof. Note first that the nonlinear Feynman Kac formula tells us that Y t = u(t,X t ) and 
Z/ = L u(t,Xt), I = 1,.. . , d, where the existence of u, Vu is guaranteed by Theorem |2.3| . 
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The proof is a simple consequence of the expansion {2.3} , using the hierarchical set ^4(3), 
applied to the integrand of the integral 



E, 



£ f(s,x s )d S y- 



f ( L a f(t i ,x i )j a m tliU + J a [L a H--x-)ku) 



du 
(3.4) 



and to f ( f; + i , X, + i 



2 



(f(x h v;-,z I ) + E i [^(x i+1 ,K + i,z ;+1 )]) 



(3.5) 



If a S ^4o(3), then one of the following holds 

(i), ( = d 

(0,0,(1,0), i = 1 c / 

{i,j),(i,j,k) i,j,k = ^ d, 

Elementarg facts about stochastic integrals and Lemma IZ2l tell us that E s [7 a [1] s u ] = 0, 
u > s for all a's as above except when a = (0), or j = 1, ... , d. For these two cases, 
we have 



E. 



[/ 0) [1 k u] =u-s and E s [7^[1 ] s , u ] = ^E [ ( W' u - W' s f ] = - 



Substituting the above into j3.4| . {3.5} and taking their difference we have 

<5,+i 



E, 



f(s,X s )ds - ( f(t h X ti ) + f[t M ,X tMi 









< 


E, 


r ( 







^ 7° [/.«?(•, X.)] tf , I du-%1 ^ J c [L°f(.,X.)] ti 



7G.Ao(5)yio(3) 



aeA(5)\A(3) 



The estimates on the remainder process complete the proof. 



□ 



The above result combined with a Stratonovich Taylor expansion applied on Z (as the 
integrand of the stochastic part) give us a first intuitive approach for the second order 
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discretization: 



Yt t * Y tM + ( f(X t „ Y t „Z tl ) + f(X tM , Y tM , Z tM ) ) 

d ( d fti+i 

- Y{ L ' U ( t h X tl ) AM// +1 + Yi L {kJ) u ( t h X ti ) / J^] ti , 5 clWi 

+ i [k ' hl) u ( t h X tl ) f ' + J {kJ) [-\}t h sdWl (3.6) 
fc,y=i Jt < 

+ /_<°'"i/(f i ,X f ,) f' + \dW l s + Y_ L a * {l] u(t h X tt ) f f ' +1 y°[1] ti . s £/W s ' 

Jti ll«|[=3 Jf ' 

+ R 3 (/-'u,^,s)c/M/ s '} 

As in (33), we need to innovate a way to recover Z/, / = 1, ... ,d, i = 0, . . . , n — 1, from 
|3.6| , but this time up to a second order error. 



If we multiply both sides of {3.6} by -g-y 1 and condition with respect to T u we shall obtain 
Z t but some surviviny terms of order 5, + i will render the approximation first order. For 
example consider 



1 



5, 



-E 



+1 



L^'hi ( t h X tl ) J s c/W s 'AW// +1 \T ti 



^L^u(t b X tl ). 



Hence, we need to find an appropriate weight, that will multiply fl3.6| and after conditioning 
with respect to T ti , provides a second order approximation for Z f . by canceling out all first 
order terms. We make the foLLowiny judicious choice 



5t 



+1 



/ = 1 d, i = n - 1 . 



(3.7) 



With a few straiyhtforwarcl computations, we obtain for any q = 1, d 

d 



E 



E 



E 



^/.'L/(t„X f ,)AM// +1 Z, q 



L /=1 

d 



Ft, 



^~Z.<°''>u(t;,X tf ) / sc/H/ s 'z; 

(=1 ' h 

Y Y ^u(t h x tl ) r ,+1 y(^[i] fi , sC /vv s 'z I (f 

f ' / d rti+1 /-f/+i 

t L / E / 



/=1 \ fc=1 || n || =3 " 

By choosing Ai = 4, A2 = —6, we have the following 



= 0. 
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Lemma 3.2. Let assumptions (A), (13) hold true and let u(t,x) denote the classical solution 
of PDE ED. Set 



,d, i = n-2. 



5, 



+1 



Then 



4 - 



<5} + , max ||/.« u (f i+1> .)|| c 

||q||— 4,5 



Proof. It should be clear from the discussion preceding Lemma [3~2l co m bine d with the result 
of Lemma [3TTI that 



5, 



< 



E, 



^ +1 R 3 (l ti.s) ds + £ +1 R 3 (L'u, tt, s)clW ; 



The result then follows from the estimate |2.5| on the remainder process and the Cauchy 
Schwartz Inequality. □ 



Le m m a s 1 3 . 11 a n cl U2\ dictate the following second order algorithm 



•Initialization 
If (C1) is In force : 

Yi,t B ■= 4>(Xn), 
If (C2) Is In force: 



:= 0, and := , := 



Z 2 \ :=V0(X n )V(X n ) 



•Backward induction: 



4!t, = Ef 



5i 



Yz^+^HX.+ u)!, ,Z] t ))Z, 
5, 



, z r .= {z) z\ 



d\T 



In = E, [ V£ ti+1 ] + J±±(f{X h Yl t ,Zl ti ) + E ( [f (X i+1 , Kf t(+1 , Z£ f i+1 )]) 



(3.8) 



A small clarification Is perhaps In order for the peculiarity of our scheme In the first back- 
ward step. If the terminal condition Is not smooth, then none of the intuitive arguments 
that we have presented can apply at time t n = T. Hence, when working under (C1), we 
take the first backward step in an Euler fashion. After this, we expect the PDE to "smooth 
out" the value function and hence all that we have discussed so far apply. We will see in 
Corollary I3.4I how a non equidistant partition can compensate for this. 
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Theorem 3.3. Let assumptions (A), (13) and either of (C1), (C2) hold true. Then there exists 
a constant C > such that 



max 

0<i<n-2 



Y tt -Y£ ti \ 2 + ^\Z tt -Z£ ti 



< c [\ - y 2 %_, | 2 + ^ | z tR _, - z 2 \_ t \ 2 )+'t_ 5 ? ™% ^ a u(t h -ni 

Proof. In the following proof, C will denote a constant whose value might change from line 
to line. It will however be Independent of the partition and of the bounds of the derivatives 
of the solution of j1 .3} . For ease of notation, we set 

Af Y := Y t , - Y 2 \ , Af Z := Z t , - Zf f| , 
Af f := f ( X t „ Y t „Z ti )-f[X t „ V£ fi ,Z£ t( ) 

:= V, +1 + ^f(X tM , Y tm , Z (i+1 ), Mif +1 := K£ tj+1 + ^f(X tM , V£ ti+1 , Z£ tj+1 ) 
AV i+1 =MJ i+1 -Mif +1 . 

Let us fix a value for i = 0, . . . , n — 2. We consider the difference of the solution of the 
BSDE at time f, and of scheme j3.8| : 



Af Y = Ei [Af +1 Y\ + %l£, [f (X h Yl h , ) + f (X i+1 , K£ ti+1 , Z&J] 



x rf, + , (3-9) 

■^-E, [? (f tl X,) + ? (t,+i , X f+ i)j - Jf 1/ [f (s, Xs)] els 



(3.10) 



According to the estimates of Lemma I37TI we have that 

Si, 



Ei [f (X b Y :,,/:,) + f VX, +1 ,Z£, + J] - -^-Ei [f(t h Xi) + f(t M ,X i+ i)] 



<C5f + , max \\L a u(t M , ^ 

||a|[— 4,5 

Moreover, according to the mean value theorem, we can argue on the existence of a real 
number and vector u<\ E R, v-\ e M. d bounded by K, such that 

^l A f/ = ^l( w Afr + v 1 -AfZ) (3.11) 
Combining |3.9 | - {3.1 1| with Young's Inequality with yi > 0, we have 

|AfK| 2 <(1+y 1 5 i+1 )|E i [Af +1 ^]| 2 + ) C5 li ( I Af V | 2 + | Af Z | 2 ) 



1 

yi<5, +1 I ~~' +[ 101=4,5' 



(3.12) 



+ ( 1 +T7^ I C5 w .maxj|/. n u(f l+1 ,-)ll^ 



Next, observe that for any random variable F which is measurable with respect to Tt M we 
have 

|E,[F2j]| 2 = \Ei[(F-Ei[F])Zi}\ 2 < J- (e, [F 2 ]-E,[F] 2 ) 
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Combining this with definition of Zf t , i = 0, . . . , n — 2 and the conclusion of Lemma [3721 we 
have 



5 /+ iE [ | Af Z | 2 ] < 2d ( E [ | A% +1 | 2 ] - E [ | E, [ AY i+1 ] | 2 ] ) + Cdf + , max. || Z.°uft + i , 



(3.13) 



Putting together (3.1 2} and |3.13| we get 

E 

< (1 +y 1 <5 i+1 )E|E i [Af +1 



|Afy|2+ ^5 L|Afz|2 



+ C«5 i+1 E[|Afy| 2 ] + |^. + ± + C5 i+1 ) 5 i+1 E[|AfZ| 2 ] 
+ C3f +1 max ||L a u(t, + i,.)||L 



(3.14) 



||a||=4,5 

< (1 + yi 5 l+1 )E[|ArV/l J | 2 ] 
C5 i+1 E[|Ary| 2 l +Caf +1 max ||£ ff o(t, + i,-)|lL 

L J || cr || =4,5 

where we have chosen yi = C4d. 

We can argue once more with the mean value theorem and Young's inequality to deduce 
that 



(1 -C5 l+1 )E 

< (1 + Cff i+ i 



Af Y\ 2 + |AfZ' 2 



4d 
I Af +1 r 



2 + ^|A^! 2 



+ C5f +1 max ||£ a u(t l+1 , 

||ar||=4 r 5 



for some different constant C. By appealing to the discrete version of Cronwall's lemma 
we complete the proof. □ 

Theorem I3.3I justifies the fact that scheme {3.8} produces a second order discretization of 
the backward component of the BSDE when (A), (B) hold true and <t> is smooth. Clearly in 
this case, max||„|i = 45 ||/_ a u(f,- + i , -JH^ is uniformly bounded, the error on the first backward 
step is of 0(1 /n ) when n points are used on the time discretization and hence, with a 
uniform partition n points, we obtain an overall error of 0(1 /n 2 ). 

The more Interesting case ocurs when <t> is only Lipschitz continuous. Given the results 
of Theorem I2.3I we expect the bound on the derivatives of u to explode as t f T. To 
compensate for this and for the fact that our first backward step is of Euler style, we work 
with a non equidistant partition that becomes more dense as we approach T. Hence we 
are still able to achieve a 1/n 2 rate of convergence with a n-points partition. 

Corollary 3.4. Let assumptions (A), (B) hold true and assume that <t> is Lipschitz continuous. 
Consider the discretization |3.8| along the partition tt : 



T 1 



1 



i = 0, 



iS > 5. 
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Then, there exists a constant C independent of the partition and of the value function u, 
such that 

|t 1 1/2 C 

max E ' 



0<i<n-1 



\Y ti -Y 2 » ti \ 2 +\Z t: -Z 2 " ti \ 2 ] <- 7 . 



Proof. Clearly, we can establish the result by estimating the terms <5 ( 5 +1 max|| a || = 4 5 \\L a u(ti + i , ■ 
According to Theorem 12.31 we have 



C||V0|| c 



1 , 5 C||V*|| c 



^5? max \\L*u(t l ,.)\\ 2 0O ^l% {J . )t 

'1-1 / /-1- — 



"- 1 T o5 d _H] 5 M 



since /3 > 5. 



To complete our proof we estimate the error on the first backward step. Using the standard 
estimate on the Euler discretization error of Lebesque Integrals we have 



E, 



[ f(s,X s )ds 
Jt„-i 



+ f ( X„_i , V„-1 , 2)7-1 ) 5„ 

2 



-f ( X n _i , , Z£ tn _ : ) <5„ | 

< cs 2 n ( 1 + 1 y n _! - I 2 + 1 z,,^ - z^, I 2 ) 



where once again, we have used the mean value theorem. Rearranging terms, we may 
argue on the existence of a constant C such that 



( 1 - C5 <> ) ( I v n -i - I 2 + % I z n -i - z^, i 2 ) 

<C5 2 n + C5 n \Z n ^-Z,\_ : | 2 
From standard estimates on BSDEs we know that under (A), (B) and (C1) 

sup E | Z ( | 2 < +00. 



(3.16) 



0<t<T 

Lastly, under (C1) we have that 



E 



[l z u,-, | 2 ] = T2 E [ I E "-i I *(X n )AVI/ n ] | 2 j 

; . . (3.17) 

= ^E [ I E„_i [ ( <D(X n ) - 4>(X n _i) ) AW„ }\ 2 \<C 
Substituting |3.17| Into |3.16| and then taking square roots, completes the proof. □ 
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